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We study spontaneously generated entanglement (SGE) between two identical 
multilevel atoms in free space via vacuum-induced radiative coupling. We show that 
the SGE in two- atom systems may initially increase with time but eventually vanishes 
in the time scale determined by the excited state lifetime and radiative coupling 
strength between the two atoms. We demonstrate that a steady-state SGE can be 
established by incoherently pumping the excited states of the two-atom system. We 
have shown that an appropriate rate of incoherent pump can help producing optimal 
SGE. The multilevel systems offer us more chanel to establish entanglement. The 
system under consideration could be realized in a tight trap or atoms/ions doped in 
a solid substrate. 

PACS numbers: 03.65.Ud, 42.50.Dv 

I. INTRODUCTION 



The recent development of quantum technologies strives to resolve the quest for the 
best entanglement source in quantum optical systems. Though entanglement is observed 
in a variety of systems, entanglement in atomic systems are favored as more scalable and 
practical systems, compared to their "photonic only" counterparts due to the development 
of reliable state-of-the-art technologies to control the atoms one-at-a-time [1] which can be 
precisely scaled to many atoms [2, 3]. Many exciting developments of entanglement sources 
are involved in atomic systems such as entanglement via atom-cavity coupling [4], atom-atom 
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entanglement via cavity [5], entanglement in trapped ions/atoms [3, 6], atomic entanglement 
in an optical lattice by controlled collision [7] and also atomic entanglement via external 
fields [8, 9]. Recently, it has also been shown that well separated atomic ensembles can 
also be entangled via coherent coupling between them [10]. Scully has extensively discussed 
entanglement in two, three and many atoms via a single photon [11] and has shown that 
such an ensemble can produce directional spontaneous emission [12]. 

Amongst the different atomic entanglement generation processes, an interesting and use- 
ful category is spontaneously generated entanglement (SGE) sources via interaction of atoms 
with a common bath of cavity field [5], vacuum [3, 6, 13], heat bath [14] or even a spin chain 
[15]. Usually, the baths have very short correlation time and hence potentially cause dis- 
entanglement [17] and decoherence [16] in an entangled system. Agarwal and Patnaik [19] 
have shown that coherences in two multilevel atoms can be generated from their interaction 
with a common vacuum bath via the retarded dipole-dipole (dd) coupling when they are 
placed in the close proximity of each other. Effect of such vacuum induced coherence (VIC) 
on the collective resonance fluorescence is discussed in [20]. SGE is particularly interesting 
from the application point of view because: practical quantum devices are often unavoidably 
coupled to the environmental bath and hence SGE can occur naturally. 

Most of the works listed above are focused on SGE in two-level atoms. Study of multi- 
level systems are important because in certain situations, participation of additional internal 
atomic levels in the process of generating entanglement or causing disentanglement is un- 
avoidable. For example, when two atoms, having a triplet P-state as their excited [19] or 
ground [20] state, are placed in close proximity, i.e., the interatomic distance is less than the 
wavelength of the atomic transitions involved, R < Xq, even the dipoles involving the or- 
thogonal o~± transitions can radiatively couple to generate additional coherences. Recently, 
Keitel and coworkers [21] have explicitly shown that the two-level approximation fails in such 
a situation. Furthermore, multilevel systems can open up new channels in bath assisted SGE 
in a vary natural way and even can give more control parameters [18]. To the best of our 
knowledge, only a few studies addresses entanglement in three-level atoms interacting with 
a continuum via the radiative coupling [17, 18]. 

In this paper, we investigate the steady state SGE between two radiatively coupled 
and incoherently pumped atoms having their energy levels in a V configuration with non- 
degenerate excited states; see Fig. 1. We derive a master equation and trace over the field 
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part to obtain the equations for the atomic dynamics. We obtain an analytical solution to 
show that multilevel systems are preferable compared to two-level systems for SGE because 
they add additional coupling channels to enhance the entanglement. We demonstrate that 
the SGE can be sustained to achieve a steady state entanglement by incoherently pumping 
the atoms. It may be noted that we are working outside the regime of VIC. In our two-atom 
system VIC could be generated if the excited states of both the atoms are degenerate or 
near degenerate [19, 20]. 

The organization of the paper is the following: In Sec. II we derive a master equation for 
the two atoms interacting with common vacuum. Tracing over the vacuum bath parameters, 
we obtain the equations for the system dynamics. In Sec. Ill, we present the time evolution 
of the entanglement between the two atoms that occurs only for a short period of time. 
In Sec. IV, we derive the atomic density matrix equations when atoms are incoherently 
pumped. We show that a steady state entanglement can be obtained between the two 
atoms purely via incoherent processes. We summarize and discuss our result in Sec. V. 




atom A atom B 

FIG. 1: (Color online) The two identical atoms under consideration. The V type atoms with non- 
degenerate excited states. The distance between the two atoms is considered to be small compared 
to the radiation wavelength, R < Xq. 
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II. THE TWO-ATOM SYSTEM AND THEIR DYNAMICS 

We consider two identical three level V systems (say A and B) in free space having two 
excited states \e a ) and \fi a ) (a = A, B), and a ground state \g a ), as depicted in Fig. 1. Both 
the atoms couple to the same vacuum field. While we do not wish to loose the generality of 
our results, our scheme can correspond to two 40 Ca atoms in a magneto-optical trap (MOT) 
in presence of a static magnetic field. The ground state can correspond to 4 1 < 5 , state and 
the excited states can correspond to the magnetic sublevels A 1 P 1 of Ca atom. The static 
magnetic field would remove the degeneracy of the A 1 P 1 sublevels and the states \e a ) and 
\ji a ) can correspond to m\ = ±1 levels. Note that we restrict to the situation where the cross 
couplings between \e a ) \g a ) and \fj, a ) <-> |e Q ) transitions are eleminated by considering 
the non-degenerate excited states. Thus the photon emitted from \e a ) — > \g a ) (|aO \da)) 
can only be absorbed by \g a ) — > \e a ) (\g a ) — > I/O)- For simplicity we consider only the case 
of real dipole moments for our discussion below. These results can be easily generalized to 
complex dipoles, e.g., involving magnetic sublevels. 

In this section, we derive the system dynamics with only the contributions from the two- 
atom coupling with vacuum, i.e., in absence of incoherent pumping. The role of incoherent 
pumping will be discussed in detail in the Sec. IV. The Hamiltonian of the two-atom system 
interacting with the vacuum field can be written in the interaction picture as 



ot=A,B ks 

+ ^2 ff akve <[p "* a+(w2 ~ w * v) * ] + ff.c.) , 

k's' 

where, the vacuum Rabi coupling coefficients corresponding to atom a are 

1/2 



(1) 



/ 2nHuj ks \ 

V v ) 



d eg = H ~^r~ ^9 ■ ?ks 



and dl g = i[ — — p^ g ■ 4v, (2) 



and (p^g) i s ^ ne dipole matrix element corresponding to the transition operator a" g = 
\ e a)(ga\ (&ug = lA t a)(fl'a|)) ^ksi^k's 1 ) is the unit polarization vector of the vacuum mode with 
frequency uiks {^k's')i a ks (a-k's') is the photon annihilation operator corresponding to the 
vacuum field with wave vector k (k f ) and polarization s (s 1 ), u>i (o; 2 ) is the atomic frequency 
corresponding to \e a ) <-> \g a ) (\fJ> a ) <-> \g a )) transitions, and x a is the position of the atom 
a. 
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We use the Zwanzig projection operator method [22, 23] to trace over the field degrees of 
freedom and obtain a reduced density matrix equation for the atoms. We use the Born and 
Markoff approximation to obtain a memoryless master equation. Referring to Ref. [19] and 
without duplicating the lengthy calculation, we write the reduced density matrix equation 
for the atoms as 

P = -i[V d d, p] + (£ s + £dd)p, (3) 

where 

Vdd = G x a% ® af e + G 2 a* ® <r* + H.c. , (4) 

is the part of the dd-interaction that contributes to the level shift. The coupling coefficients 
are 



ks 



(6) 



Here R = xa — %b- Further, the Liouvillian operators £j are: 

£ s p = 7i[(2< e p< - af eP - pafj + A —> B] 
+ 7 2 [(2<P< - <r> - P<J + A -> B], 
corresponding to spontaneous emission of the atoms, and 

CddP = [Vi{2a B ye pai g - af g af e p - pa%a%) + H.c] 

+ [T 2 (2af^ g - a^p - p«) + H.c], 

corresponding to the dd coupling mediated by the vacuum. Note that the subscript in p a is 
dropped for brevity. Here the spontaneous decay rates are given as 



(7) 



ql 

and the atom-atom coupling coefficients are obtained as 



7 1 = /^S 71 "^ 1 -^fcOKJ 2 , 

ks 

7 2 = ^^^(^-^JqdKgl 2 , (8) 



ks 

r 2 = ^E^-^ld^re* 5 - (9) 
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Further the index a has been dropped as we consider that the atomic dipoles corresponding 
to the same atomic transitions are parallel to each other, i.e., pfj || pfj. Clearly, the radiative 
coupling terms Tj and G{ have numerical significance only in the limit \kR\ and \k'R\ < 1. In 
the other limit, when the interatomic distance R is too large, only the spontaneous emission 
terms survive and the atoms behave as two independent atoms. We refer to [19] for the 
detail steps of the calculation. 

Now let us assume that initially, the two-atom state is |ep), where the notation \ij) = \i A )® 
\3b), i,j = e, fj,,g. The nine two-atom basis states are |ee), \e/j), \eg), \pe), \pp), \pg), \ge), 
\gp) and \gg), and we number them 1 through 9 in the same order as above to simplify the 
notations for the density matrix elements (^aJbIpI^l/b)- For examples, the density matrix 
element (eAPB\p\ e APB) corresponding to our initial state | e/x) is represented as p 22 in the 
new notation. The full density matrix equation involves 81 matrix elements but for the 
above initial condition, many elements would be identically zero and only 10 density matrix 
elements survive. We consider the geometry where dipole matrix elements are orthogonal to 
each other and are real (as discussed in Sec. Ill of [19]), such that the parameters G i} 7 i5 
are real numbers seen Eq. (26) in [19]. In the following, we explicitely write the dynamics 
equations only for those surviving density matrix elements as 

P22 = -2(7i + 7 2 )p 22 , 

P33 = - 2 7lP 3 3 + 2 7 2 P22- r i(P73 + P37) 

-iGi(p 73 - p 37 ), 



P37 = 


- 2 7lP37 


- Ti(p 77 + P33) - 


iGi(p 77 


-P3s)> 


P66 = 


~ 2 7 2 P 6 6 


- r 2 (p 86 + p 68 ) - 


iG 2 {p m 


-P 6 s) ? 


Pes = 


- 2 7 2 P68 


- r 2 (p 88 + p 66 ) - 


iG 2 {p 8 8 


-Pee). 


P77 = 


-2 7l p 77 


- r i(p 3 7 + P73) - 


iGi(p 37 


-P73)> 


P88 = 


- 2 7 2 Pss 


+ 2 7ip 22 -r 2 (p 86 


+ P 6 s) 





-iG 2 (p 68 - p 86 ), 

P99 = 2 7l(P33 + P77)+ 2 7 2 (p66+P88) 

+2ri(p 37 + p 73 ) + 2r 2 (p 68 + p 86 ). 

Note that the conjugate matrix elements p 73 and p 86 (conjugates of p 37 and p 68 , respectively) 
also evolve. Using the Laplace transform method, we solve the above coupled equations for 
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the density matrix elements with the initial condition p 22 = 1 to obtain their time evolution 



p 22 (t) = e- 2 ^)', 

"7 2 cosh(2rit) - r\ sinh(2rit) - -f 2 e' 2 ^ t 
7 2 cos(2G 1 t) + G x sin(2G^) - 7 2 e" 27 2*- 

+ iTTgi J' 

P37W = x 

-r ie - 2 T2* _ r x cosh(2rit) + 72 sinh(2rii) 

[ fF^I 

iG x cos(2G^) - 2 72 sin(2G 1 t) - zde" 2 ^*! 

GfT^l J' 

Tie" 272 ' 
PeeW = iL ^— x 

"7! cosh(2r 2 t) - T 2 sinh(2r 2 t) - 7 1 e- 27 i* 

L /l 1 2 

7i cos(2G 2 t) + G 2 sin(2G 2 t) - 7ie~ 27lf 1 
7? + G 2 J 
Pea(t) = ptAt)]^, p 77 (t) =p 66 (*)| 1 „ 2 , 

P 88 W = P33(*)|i„ 2 ' 

P99W = 1 -P33W -Pee(*) -P77W -Pk»(*)- 

It may be noted that the initial state p 22 = (eAP_B|p|eyi^ B ) decays with a rate of the sum 
of the decays of both excited states but does not depend on the dd coupling terms Tj and 
Gj. However, the other population and cross terms strongly depend on the dd coupling. 
The dd terms Tj play the role of decays via the cosine and sine hyperbolic functions and the 
coefficients Gi cause the dd coupling induced vacuum Rabi oscillations. The time dependent 
solutions of the matrix elements show the oscillations with frequencies determined by the 
atom- atom coupling coefficients G\ and G 2 - Further, the dd terms are strongly dependent 
on the interatomic distance R. Hence the dynamics of the density matrix elements are also 
strongly affected by R. We will present numerical plots and discussions for some of the 
important density matrix elements that help in evolving the SGE in the following section. 
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III. TIME EVOLUTION OF SGE 

In this section, we calculate the time evolution of the entanglement between two-atoms. 
Out of various different methods to calculate entanglement between the two atoms, we 
choose the negativity, defined by [24] 

N(P)= hTA \~ l = Y1 ' A - ( 12 ) 

as the measure of entanglement. Here ||p TA || denotes the trace norm of p TA [25]; p TA is 
the partial transposition matrix of the atomic system density operator p{t). The primed 
sum in the above equation represents the sum over only the negative eigenvalues Aj of p TA . 
For a high dimensional system, while a non-zero N(p) is a sufficient condition to prove 
that a system entangled, but a null N(p) does not necessarily qualify a system to become 
disentangled. From the definition, the negativity N can also be greater than 1. For different 
dimensions of the density matrix, the maximum value of N is different. For a two- atom 
three-level system, such as ours, the state \& = ^[|ee) + \pp) + \gg) is maximally entangled 
one with the negativity N(^>) = 1. 

Thus to obtain negativity in our two-atom system, we first calculate the eigenvalues of 
p TA by using the density matrix p(t). After a lengthy calculation, we obtain the exact 
eigenvalues as 

Ai = A 2 = 0, A 3 = p 22 (t), A 4 = p 33 (t), 
A5 = Pm(t), A 6 = p 77 {t),X 7 = PssW, 

As = ^ + + 4[|p 37 WI 2 + Mm, (13) 

A 9 = ^-y P Ut)+A[\p 37 (tW + \p 68 m- 

It is clear that among all of the above eigenvalues, only A 9 can become negative. Therefore, 
if we obtain a negative A 9 , it is sufficient to prove the occurence of the SGE. Clearly, A 9 can 
be negetive only if at least one of the matrix element p 37 = (eg\p\ge) or p 68 = (pg\p\gp) is 
non-zero, i.e., if there is an exchange of at least one photon between two atoms. Thus it is 
clearly established that radiative coupling leads to SGE in the two-atom system. Because 
both p 37 and p 68 contribute to the generation of entanglement. Thus the three-level atoms 
offer us more channels to establish entanglement than two-level one. 
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FIG. 2: (Color online) The time-dependence of real and imaginary parts of matrix element 
/9 37 (t) with different values of Ti(r2) and G\(G2)- We set 7 = l,r = 1.2. The parameters 
= (2.4,0.9), (0.9,0.8), and (-0.24,0.2) correspond to R = 0.83Ai, I.I8A1 and 2.78Ai that 
are represented in the figure as the solid, dashed and dotted line, respectively. 

Before we proceed further, we first study the density matrix elements p 37 (t) and p 99 (t) 
that determine the entanglement. Note that the two excited states are non-degenerate. 
Assuming that cu 2 = ru>i, and say j 1 — 7; Ti — T, G\ = G, we have 7 2 = rj, T 2 = rT 
and G2 = rG. For the numerical plots presented below, we use the parameters given in the 
Table I determined from their definitions in Eqs. (5, 8, 9) [30]. Note that as the interatomic 
distance reduces, V oscillates [19]. However, peak value of V decreases with larger interatomic 
distances due to decreased dipole-dipole coupling between the two atoms. Thus, Tj ( Gi) 
can be same in two or more interatomic different distances. But any particular interatomic 
distance determines a particular value for the pair of (r^, Gi). We will concentrate on the the 
photon exchange process with decreased trend of Tj ( Gi ) with the increasing interatomic 
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separation. So, the parameters given in table is approximately monotonic. 

TABLE I: The interatomic distance and the corresponding coupling parameters. All the frequency 
units are scaled with 7. 



R (in unit of Ao) 


Ti 


r 2 


Gi 


G 2 


0.50 


0.967! 


O.9672 


8.07! 


8.0 72 


0.83 


0.97i 


0.9 72 


2.4 7 i 


2.4 72 


1.18 


O.871 


0.8 72 


O.971 


0.9 72 


2.78 


0.2 7 i 


0.2 7 2 


-0.247! 


-0.24 72 



In Fig. 2, we present the evolution of the density matrix element p 37 = (eg\p\ge) repre- 
senting the single photon radiative coupling process, in which atom A loses its excitation to 
excite atom B from its ground state \g&) to the state \e B )- It is observed that, initially the 
process of exchange of photon increases. However, after reaching a maximum, Rep 37 falls 
off quickly within one spontaneous emission cycle. The time needed to reach the maximum 
value for Rep 37 is determined by r -1 . The maximum of Imp 37 occurs at t ~ 7 _1 . In long 
time limit both real and imaginary part of p 37 vanish. Similar conclusions can be derived 
for the matrix element p 68 which physically represents the simultaneous probability of two 
processes \p A ) -> \g A ) and \g B ) -> \p B )- 

In Fig. 3 we present the evolution of the population in the state \gg), i.e., the matrix 
element p 99 (t). This plot also supports the physical process we described above. We can 
see that again large values of T slows down p 99 {t) to reach at its asymptotic value. In other 
words, the radiative coupling process survives longer. The steady state value is p 99 = 1, i.e., 
both the atoms reach their ground states in the long time limit. 

Next, we discuss the property of the entanglement generated this atomic system. We 
present the plot of N(t) that describes the time evolution of the SGE for different values 
of T and G in Fig. 4. It is shown that at t — 0, there is no entanglement because initially 
the atomic system is in the state \ep). For t > the entanglement evolves to reach its 
maximum value, and then undergoes a process of disentanglement. Finally, steady state 
of negativity becomes identically zero. But the relaxation time becomes longer when the 
atoms are nearer (seen solid line in Fig. 2a and Fig. 4 ). From the solution Eq.(ll), 
we know that relaxation time of Re (p 37 ) is determined by max{^^-^, 2{^~r\)} which 
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2 4 6 8 10 12 

FIG. 3: (Color online) The matrix element for our atomic system as a function of t corresponding 
to different pairs of (G, T). The curves and also corresponding legends are same as in Fig. 2. 

means larger the value of longer the relaxation time for disentanglement. This is the 
competition between dd coupling and the all-direction spontaneous emission. The evolution 
of SGE presented here can be understood as due to the following multiple radiative exchange 
processes: spontaneous emission of atom A via |e^) — > (atom B via — > |<?b)) 
transition is followed by an absorption of the spontaneously emitted radiation by atom B 
via \qb) — > \gb) (atom A via \qa) — > l/^)), as it is expressed by p 37 (£) (P6sW) causing the 
two- atom entanglement. 
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FIG. 4: (Color online) The negativity for our atomic system as a function of t corresponding to 
different pairs of (G, T). All the parameters correspond to Fig. 2, respectively. 

IV. STEADY STATE ENTANGLEMENT WITH BROADBAND INCOHERENT 

PUMPING 



From the previous section, we have seen that the disadvantage of SGE is the temporal 
evolution and the quickly diminishing of SGE due to the decays in the system. However, 
for any useful application, for example, to use the two-atom system as a coupled qubit, 
sustaining the generated SGE is essential. In order to achieve a steady state entanglement, 
we introduce an incoherent pump to continually repump population to the excited states. 
Steady state entanglement using a classical coherent pumping is discussed in [26] where the 
coherent pumping with single frequency drive the atoms. Coherent field assisted entangle- 
ment is rather intutive. On the contrary, the usual notion associated with an incoherent 
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pumping with a decoherence process and hence a source of disentanglement. However, in 
what follows below, we show that an appropriate strength of incoherent pump can lead to a 
steady entanglement between the two atoms. 

We consider a broadband incoherent pump acting on the two atoms which incoherently 
drive the population from \g) to |e) (|/i)) levels. Incoherent pumping can be modeled as an 
inverse process of spontaneous emission [27, 28]. Thus we add a third Liouvillian £ inc to the 
the master equation Eq.(4) to get 

P = -i[Vdd, p] + (£ s + Cm + AncV, (14) 

where 

C mcP = A 1 [(2af gP af e - a A gg p - po A gg ) + A -> B] 

+A 2 [(2aX - a^p - pa A J + A - B\. 

(15) 

where Ai and A 2 denote incoherent pumping rates for \g a ) — > \e a ) and \g a ) — > \p a ) transitions 
respectively. We explicitly write the density matrix equations involved in presence of the 
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incoherent pump as 

Pn = - 4 7iPn + 2Ai(p 77 + p 33 ) 

P22 = -2(7i + 7 2 )p 22 + 2Aip 88 + 2A 2 p 33 , 

P33 = -2siP 3 3 + 272P22 + 271P11 + 2Aip 99 

-« G i(p73 - P37) - r i(p73 + P37)) 
P37 = -2sip 37 - Ti(p 77 + p 33 ) 

+2rip n -«Gi(p 77 -p 33 ), 
p 4 4 = -2(7i+7 2 )p44 + 2Aip 66 + 2A 2 p 77 

P55 = -47 2 P 5 5 + 2A 2 (p 66 + p 88 ) (16) 
Pee = -2s 2 p 6 6 + 2 7lP44 + 272P55 + 2A2P99 

-iG 2 {p m - p 68 ) - r 2 (p 86 + peg), 
Pes = -2s 2 p 68 - r 2(p 88 + Pee) + 2rip 55 - iG 2 (p 88 - p 66 ), 
p 77 = -2 Sl p 77 + 2 7l p n + 2 72 p 44 + 2A!p 99 

-r\(p 37 + p 73 ) - iGi(p 37 - p 73 ), 

P 88 = -2s 2 p 88 + 2 7l p 22 + 2 72 p 55 + 2A 2 p 99 

-r 2 (p 86 + p 68 ) - iG 2 (p 68 - Pse)) 
P99 = 2 7l (p 33 + p 77 ) + 2 72 (p 66 + p 88 ) - 4(Ai + A 2 )p 99 

+2r 1 (p 37 + p 73 ) + 2r 2 (p 68 + p 86 ). 

where s a = 7a + A 1 + A 2 (a = 1,2). Note that in presence of the incoherent pumping Aj, 
the non-zero additional terms are the populations \ee), |pp), |pe). Hence we have a total 
of 13 non-vanishing density matrix elements in presence of Aj. Since we are looking for a 
steady state SGE, we calculate the steady state values of the matrix elements by setting the 
differentials in the left-hand sides as zero and solving the coupled equations. The analytical 
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solutions that we obtain are 



Pn = y A i a 2, 

P22 = 77 — \ 2 v (aiAi + a 2 A 2 ), 
&(7i + 7 2 ) 
7i7 2 

P33 - ~~ ^ ° 2 ' 

P37 = ( Al ~~ 7l)«2, P44 = P22, 

7iOi a _ 7i7 2 a i 

P55 - ~~g~ A2 ' ^66 - ^ ' 

T1 ^2 

Pes = ^-(A 2 -72) a i 5 

P77 = P33> P88 = P66; 

P99 = ^f 1 (f3 1 a 2 + a 1 (3 2 ), (17) 



with 



ai = A 2 [7i + 2/3i(7i + 7 2 )], 
a 2 = A 1 [ 72 + 2/3 2 ( 7l + 7 2 )], 

„ ^7' + r^(A a - tJ 

^ 2 Sa7a (A 1 + A 2 ) ' J 
6 = 2 7l72 (7i + I2W1 + 1)«2 + (/3 2 + l)ai] 

+ (7i + 72)(a272 A i + a i7l A 2 ) 

+7 1 7 2 (2a 2 A 2 + 2aiAi) (18) 

It is interesting to note that the steady state of the density matrix elements do not depend 
on G\ and G 2 . The level shift parameters G\ and G 2 typically contribute to oscillation of 
population and coherence terms. Hence, in the long time limit, such fast oscillation terms 
vanish. Thus steady state solutions in Eq. (17) are independent of G\ and G 2 . 

Once again, as in the previous section, we calculate the eigenvalues of p Tjv to measure the 
entanglement. We obtain equation for the non-zero eigen values as 

(Pll - A)(p 5 5 - ^)(P99 - A) - |P 3 7| 2 (P55 ~ A) 

-|P 68 | 2 (Pii-A) = 0. (19) 

We obtain the numerical values of A solving the above equation and substitute in Eq. (12) to 
obtain the steady state negativity as a function of A;, as shown in Fig. 5. We have scaled the 
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incoherent pumping rate Aj with the spontaneous decay rate r y 1 — 7 and also for simplicity 
we have assumed A 1 = A 2 = A. Clearly, a non-zero steady state entanglement is obtained 
by incoherently repumping the excited state. Smaller the interatomic distance, larger is 
the steady state entanglement. Further, as the incoherent pumping rate is increased, the 
SGE increases but after reaching a certain optimal value at around A = O.O87, the atomic 
entanglement starts to reduce. For smaller the interatomic distances, even stronger inco- 
herent pumping can be used to get entangled atoms. Without any incoherent pumping the 
steady state SGE is identically zero. Physically, the increase in entanglement with the in- 
coherent pumping can be understood as follows: the spontaneous emission in either of the 
two atoms followed by exchange of photon between them generates SGE. But that does not 
survive longer because the spontaneously emitted photon can escape in any arbitrary direc- 
tion. Once both atoms loose their excitation, SGE vanishes. An incoherent pump assists 
the atoms to bring back to the desirable excitation so that more spontaneous emissions and 
hence photon exchanges can take place between the two atoms. Thus, increasing the re- 
pumping via incoherent pumping helps increasing the SGE. However, incoherent repumping 
also competes with the photon exchange process to re-excite the atoms. While an excitation 
due to the photon-exchange process enhances the entanglement, an excitation by incoherent 
process has no direct contribution to the entanglement. In fact, for a larger A, the incoherent 
excitation dominates the photon exchange process and hence causes a decrease in SGE. For 
A 3> T, h Gi, SGE becomes identically zero. 



V. DISCUSSION AND SUMMARY 



We have investigated the spontaneously generated entanglement in a system of two three- 
level atoms are coupled to the common vacuum field. We have presented the time evolution 
of SGE due to the photon exchange between the two atoms. We have shown that both the 
magnitude of entanglement and the survival period of SGE are enhanced by reducing the 
interatomic distance. From our analytical calculations, we have shown the strong dependence 
of the SGE on the radiative coupling parameters. We have explicitly demonstrated that the 
multilevel atoms are preferable compared to their two-level counterparts for SGE, because 
each channel adds to enhance the magnitude of the entanglement. In the long time limit, 
however, SGE vanishes. 
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0.025r 




FIG. 5: (Color online) The steady state negativity for our atomic system as a function of A, where 
we set Ai = A2 = A, 7 X = 7 = 1, r = 1.2, 7 2 = rj. From bottom to top, the values of V = 0.8, 0.9 
and 0.96, corresponding to R = I.I8A1, 0.83Ai and 0.5Ai, respectively. 

Further, to reinforce the above short term evolution of SGE in the radiatively-coupled 
two-atom system, we have proposed to use an incoherent pump that assists in repumping 
the deexcited atoms and sustain the SGE. We have demonstrated that for a certain range 
of incoherent pumping, the steady state value of SGE increases as it prevents atoms from 
loosing their excited state population. However, since incoherent pumping competes with 
the two-atom photon exchange process to reexcite the atoms, a stronger incoherent pumping 
is shown to be undesirable. We have shown that an appropriate rate of incoherent pump 
can help producing optimal SGE. 



18 



The above entanglement can further increase (not discussed here) if one considers atoms 
having degenerate or near degenerate excited states in their excited state that has additional 
coherences [19], which will be discussed elsewhere. The radiative coupling discussed above 
can be realized in any a tight ion trap. However, this work can be generalized to realizing 
SGE in a chain of quantum dots or even in a typical dense multi-atom system. We believe 
this work will open up a new way to utilize the naturally occurring SGE to realize an efficient 
entanglement source. 

This work was supported by NSFC under Grant No. 10774020, and also supported by 
SRF for ROCS, SEM. AKP is indebted to Prof. G. S. Agarwal, Dr. P. Anisomov, and Prof. 
M. O. Scully for discussions on various aspects of the two-atom multilevel systems. 



[1] B. B. Blinov, D. L. Moehring, L.-M. Duan, and C. Monroe, Nature 428, 153 (2004). 

[2] Y. Miroshnychenko, W. Alt, I. Dotsenko, L. Fster, M. Khudaverdyan, D. Meschede, D. 

Schrader, A. Rauschenbeutel, Nature 442, 151 (2006). 
[3] H. Hfner, W. Hnsel, C. F. Roos, J. Benhelm, D. Chek-al-kar, M. Chwalla, T. Krber, U. D. 

Rapol, M. Riebe, P. O. Schmidt, C. Becher, O. Ghne, W. Dr, and R. Blatt, Nature 438, 643 

(2005). 

[4] S. Haroche and J.-M. Raimond, Exploring the Quantum, Oxford Univ. Press, NY (2006); M. 

B. Plenio and S. F. Huelga, Phys. Rev. Lett. 88, 197901 (2002); M. karja, N. M. Bornik, M. 
Lffler, and H. Walther, Phys. Rev. A 60, 3229 (1999). 

[5] J. I. Cirac and P. Zoller, Phys. Rev. A 50, R2799 (1994); E. Hagley, X. Maitre, G. Nogues, 

C. Wunderlich, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 79, 1 (1997); M. 
B. Plenio, S. F. Huelga, A. Beige and P. L. Knight, Phys. Rev. A 59, 2468 (1999); A. Beige, 
S.Bose, D.Braun, S. F. Huelga, P. L. Knight, M. B. Plenio and V. Vedral, J. of Mod. Opt. 
47, 2583 (2000). 

[6] G. K. Brennen, I. H. Deutsch, and P. S. Jessen, Phys. rev. A 61, 062309 (2000). 
[7] O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. Hsch, and I. Bloch, Nature 425, 937 
(2003). 

[8] C. Cabrillo, J. I. Cirac, P. Garc-Fernn dez, and P. Zoller, Phys. Rev. A 59, 1025 - 1033 (1999); 
zgr kir, Alexander A. Klyachko, and Alexander S. Shumovsky, Phys. Rev. A 71, 034303 



19 



(2005) . 

[9] M. Kiffner, J. Evers, and C. H. Keitel, Phys. Rev. A 75, 032313 (2007). 
[10] L.-M. Duan, M. D. Lukin, J. I. Cirac, and P. Zoller, Nature (London) 414, 513 (2001); C. 

W. Chou, H. de Riedmatten, D. Felinto, S. V. Polyakov, S. J. van Enk, and H. J. Kimblel, 

Nature (London) 438, 828 (2005). 
[11] M. O. Scully, Laser Phys. 17, 635 (2007). 

[12] M.O. Scully, E. S. Fry, C. H. Raymond Ooi and K. Wodkiewicz, Phys. Rev. Lett. 96, 010501 

(2006) . 

[13] U. Akram, Z. Ficek, and S. Swain, Phys. Rev. A 62, 013413 (2000); L. Jakobezyk, J. Phys. 

A: Math. Gen. 35, 6383 (2002); Z. Ficek and R. Tanas, Phys. Rev. A 74, 024304 (2006). 
[14] M. S. Kim, J. Lee, D. Ahn, and P. L. Knight, Phys. Rev. A 65, 040101 (2002); S. Schneider 

and G. J. Milburn, Phys. Rev. A 65, 042107 (2002); J. Shang and H. Yu, Phys. Rev. A 75, 

012101 (2007). 

[15] D. Rossini, T. Calarco, V. Giovannetti, S. Montangero, and R. Fazio, Phys. Rev. A 75, 032333 

(2007) . 

[16] B. Duplantier, J.-M. Raimond, V. Rivasseau (Eds.), Quantum Decoherence, Progress in Math- 
ematical Physics, Vol. 48, Springer Verlag (2007). 
[17] L. Derkacz and L. Jakobezyk, Phys. Rev. A 61, 062309 (2000). 
[18] M. Kiffner, J. Evers, and C. H. Keitel, Phys. Rev. A 75, 032313 (2007). 
[19] G. S. Agarwal and A. K. Patnaik, Phys. Rev. A 63, 043805 (2001). 

[20] J. Evers, M. Kiffner, M. Macovei, and C. H. Keitel, Phys. Rev. A 73, 023804 (2006); M. 

Macovei, Z. Ficek, and C. H. Keitel, Phys. Rev. A 73, 063821 (2006). 
[21] M. Kiffner, J. Evers, and C. H. Keitel, Phys. Rev. A 76, 013807 (2007). 
[22] L. Mandel and E. Wolf, Optical coherence and quantum optics, Cambridge Univ. Press (1995). 
[23] G. S. Agarwal, Quantum Statistical Theories of Spontaneous Emission and Their Relation 

to Other Approaches, Spring Tracts in Modern Physics: Quantum Optics (Springer- Verlag, 

Berlin, 1974) 

[24] G. Vidal, R. F. Werner, Phys. Rev. A 65, 032314 (2002). 

[25] K. Zyczkowski, P. Horodecki, A. Sanpera and M. Lewenstein, Phys. Rev. A 58, 883 (1998). 
[26] 6. Cakir,A. A. Klyachko, and S. Shumovsky, Phys. Rev. A 71,034303(2005). 
[27] Y. Mu and C. M. Savage, Phys. Rev. A, 46, 5944 (1992). 



20 



[28] T. Salzburger and H. Ritsch, Phys. Rev. Lett. 93,063002(2004). 

[29] D. Bru/3 and C. Macchiavello, Phys. Rev. Lett. 88, 127901(2002); D. Kaszlikowski, P. 

Gnacihski, M. Zukowski, W. Miklaszewski, and A. Zeilinger, ibid. 85, 4418(2000). 
[30] The numerical parameters are obtained for a chosen geometry with 6 = ^, (p = j in Eq. (26) 

of [19]. 



